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Abstract 

A functional method is discussed, where the quantum fluctuations of a theory 
are controlled by a mass parameter and the evolution of the theory with this 
parameter is connected to its renormalization. ft is found, in the framework of 
the gradient expansion, that the coupling constant of a N = 1 Wess-Zumino 
theory in 2+1 dimensions does not get quantum corrections. 



1 Introduction 



The understanding of the notion of renormalization had done a great step for- 
ward when the connection with the blocking procedure was made [1]. In this 
framework, the control of the quantum fluctuations can be performed by the 
progressive elimination of the Fourier components of the fields in a given theory 
[2] . The parameters of the theory are then functions of the running cut-off and 
the corresponding renormalization flows describe explicitly the dependence of 
the theory on the energy scale. There are different approaches here, due to the 
freedom in defining the coarse-graining procedure. 

So as to avoid the dependence on a specific blocking procedure, an alternative 
approach to renormalization was proposed [3] . There the quantum fluctuations 



f 



are controlled by the mass of a scalar field: a large mass freezes the quantum 
corrections, the system is classical, and as the mass decreases the fluctuations 
appear and generate the quantum effects. The evolution of the parameters 
with this mass can thus be seen as renormalization flows and it was shown that 
these correspond to the usual flows at the one-loop level. Beyond one-loop, the 
"mass-controlled" flows do not coincide anymore with the loop expansion since 
they are generated in the framework of the gradient expansion, which is based 
on an expansion of the effective action (proper graphs generating functional) in 
derivatives of the field, rather than based on an expansion in powers of h. An 
example of the difference with a loop expansion is given in [4] , where the same 
method is applied to QED and the running mass is the fermions one. There the 
famous Landau pole is recovered at one-loop, but disappears if all the quantum 
corrections are taken into account. It should not be forgotten that this result 
holds at the lowest order in the gradient expansion and higher order derivatives 
were not taken into account. Finally, this functional method was also applied to 
the study of the Coulomb interaction in an electron gas [5] where the one-loop 
Lindhard function and screening are reproduced. 

In the present paper, we propose to make use of this functional method in 
2+1 dimensions, motivated by relativistic-like effective theories for the descrip- 
tion of high-temperature planar superconductors [6] . The use of renormalization 
procedures here is essential: starting from lattice models dealing with the initial 
strongly correlated electron system, the corresponding theory in the continuum 
should give to the fore the relevant operators and interactions, and the present 
method is believed to be well appropriate for this. 

Section 2 presents the method in the simple case of a self-interacting scalar 
field, so as to concentrate on the physical meaning of the flows that are obtained. 
The connection with the usual renormalization flows is made and the technical 
details, as well as the derivation of the evolution equation, are given in the 
Appendix A. 

Section 3 deals with the supersymmetric generalization of the previous case. 
N = 1 supcrsymmctry in 2+1 dimensions docs not have non-renormalization 
theorems, since the integrations over superspace have all the same measure 
d 3 xd 2 9 and no chiral supcrficld can be defined. It is thus interesting to make 
use of this functional method to learn about the behaviour of the parameters 
under the influence of quantum corrections. It will be found that, to the order 
of the gradient expansion and the truncation of the potential which are consid- 
ered, the coupling constant does not get quantum corrections, at any order in 
h. This absence of renormalization of the parameters, though, should not be 
confused with a non-renormalization theorem: higher powers of the superficld 
or higher order terms in the gradient expansion (higher derivatives, derivative 
interactions) would lead to a renormalization of the coupling constant. For the 
sake of clarity in the presentation, the derivations are given in the Appendix B. 

Finally, the conclusion contains a discussion on the gradient expansion ap- 
proximation. 
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2 Self interacting scalar field 



We present here the functional method in the case of a self interacting scalar 
field. The difference with 3+1 dimensions is that the integrals do not need 
regularization, as will be seen. 

The starting point is the Euclidean action 

S = + (1) 

where the mass dimension of the coupling is [eo] = 1 . The classical system will 
correspond to A — > oo, where the mass is infinite and thus the quantum effects 
are not present. These will gradually appear as A decreases, and when A — ► 1 
the full quantum corrections will be present. 

As usual, we introduce the effective action T as the Legendre transform 
of the connected graphs generator functional. After some manipulations (see 
Appendix A), we find then the following exact evolution equation with A: 




To obtain informations on the physical processes, we express the functional T 
in terms of its argument <j> via a gradient expansion and a truncation of the 
potential, to assume the following functional form: 

r = + + (3) 

In this approximation, we obtain the following evolution of the parameters (see 
Appendix A): 



cMAm 2 ) = mg(l 



3e 2 mg 

32^A3/2 TO 3 • W 

One should bare in mind that the equations (4) contain all the quantum 
corrections, in this approximation of the gradient expansion. The connection 
with the one-loop results is obtained by making the replacements m 2 — > m" 2 and 
e — > e in the right hand sides of Eqs.(4), since restoring the factors h leads to 
e — > he: 



d x (Xm 2 ) = ml-h^^ + 0(h 2 ) 
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The solutions are found by integrating from A — oo to A = 1: 

2 2 

m — m Q 
e-E 

where M 2 and E are constant of integration. The solutions (6) coincide with 
the usual one-loop results obtained when computing the appropriate Feynman 
graphs (see Appendix A) , where the scale M is linked to the cut-off of the theory 
and E — e . Therefore it has been shown that the evolution with the parameter 
A controls the quantum fluctuations and is consistent, at one loop, with the 
usual renormalization scheme. 

3 N = 1 Wess-Zumino model 

We deal here with the supersymmetric generalization of the previous example. 

In terms of the real scalar superfield Q, the supersymmetric generalization 
of the previous section is the Wess-Zumino model, with the following classical 
Euclidean action 

S[Q] = J d 5 z {\qD 2 Q + ^m Q 2 + |q 3 } , (7) 

where z = (x, 9) is the superspace coordinate and the conventions are those 
taken in [7]. The coupling g has mass dimension [go] = 1/2 and the on-shell 
Lagrangian contains, among other terms, the interaction g 2 (f) 4 , where <j> is the 
scalar component of Q. 

The evolution equation of the proper graphs generator functional T has a 
similar form as the one obtained in the non-supersymmetric case and is 

In the framework of the gradient expansion, and truncating the potential to the 
order Q 3 , we consider the following ansatz for the functional dependence of T: 

T[Q}= J d 5 z^QD 2 Q + ^mQ 2 + ^Q 3 Y (9) 

In this approximation, the trace of the operator (i5 2 r) _1 is computed in the 
Appendix B and it is shown that the relevant quadratic term, as far as the 
evolution of the mass is concerned, is 



= -M 2 + + o(h 2 ) 
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/d 6 k 
j^dH Q(k,6)Q(-k,0) (10) 

d 3 p [p 2 - 3(Am) 2 ] 



(2tt) 3 [p 2 + (Am) 2 ] 2 [(fc - p) 2 + (Am) 2 ] ' 

In the approximation (9), the quantum corrections to the mass are obtained 
when k = in the integrand. The identification of the non-derivative quadratic 
terms in both sides of the evolution equation (8) gives then 



d 3 p p 2 - 3(Am) 2 
(2tt) 3 [p 2 + (Am) 

The integration over p is easy and gives a vanishing integral, such that 



d x (\m)=m +g m / 7^35 ^ I ^^213 ■ ( n ) 



5 A (Am) = m , (12) 

i.e. the evolution of (Am) is classical and does not have quantum corrections. 
The evolution of m with A is then 

m = m (l + ^j , (13) 

where c is a dimensionless constant, which is computed at the one- loop level 
in the Appendix B, using the Feynman graph technique. We can note that 
in a numerical resolution of the differential equation for the mass, the large 
but finite initial value for A would not allow us to set exactly m = m since 
quantum fluctuations, even tiny, are present for finite A. This is the reason why 
the constant of integration c would not vanish. 

Let us look at the evolution of the coupling constant g(X). It is seen in the 
Appendix B that the relevant cubic term, as far as the evolution of the coupling 
is concerned, is 



g 3 J (0 (0 dH e )Q^ °)^- k <?> *) ( 14 ) 

d 3 p 4(Am)[(Am) 2 - p 2 } 

(2tt) 3 [p 2 + (Xm) 2 ] 2 [(p + q) 2 + (Am) 2 ] [(p - k) 2 + (Am) 2 ] ' 

In the framework of the approximation (9), the quantum corrections to the 
coupling are obtained when k = q = in the integrand. The identification of 
the non-derivative cubic terms in both sides of the evolution equation (8) gives 
then 



*»= i2 » 3 <W#^& (15) 
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Once again, the integration over p is easy and gives a vanishing integral, such 
that 



d\g = 0, 



(16) 



i.e. the coupling constant does not evolve with A. This is consistent with the 
one-loop result (see Appendix B) and provides a generalization to all orders in 
h, in the approximation (9). 

Finally, the quantities Am and g are frozen to their initial value given when 
A — > oo, i.e. when the quantum fluctuations can be neglected. The coupling 
constant g, in the present approximation (9), does not get any quantum correc- 
tion. 

4 Conclusion 

To conclude, let us stress again the non-perturbative nature of the flows that 
are obtained with this method. The evolution equations of the coupling con- 
stants contain all the quantum corrections, independently of any pertubativc 
expansion in K, in the approximation of the gradient expansion and the poly- 
nomial truncation of the potential. This feature is interesting in the context 
of low-energy effective theories, where only low powers of the momentum are 
kept into account in the functional dependence of the effective action, as well 
as relevant operators only, in a renormalization group sense. 

Let us be more precise concerning this gradient expansion approximation. 
The evolution equation of the effective action, giving d\T, is exact and contains 
all the quantum corrections. However, it looks like a one-loop equation, since it 
has the form (when the factor K is restored) 



where T is some operator applied on T. This is the reason why the usual one- 
loop results are obtained, when T is replaced by the classical action S in the 
right-hand side of Eq.(17). The next step forward, beyond one-loop, is to as- 
sume a functional dependence of T and to plug it inside the evolution equation 
(17). The resulting couplings' evolution consist then in a partial resumation 
of the graphs. In the present example, where no new vertices were considered 
compared to the classical ones, the resumation takes into account the one-loop- 
likc graphs, with the internal lines being the full propagators (full in the present 
approximation of the gradient expansion and polynomial truncation of the po- 
tential). 

Note that more general gradient expansions could be considered, besides 
including higher order powers of the field: one could add higher order derivatives 
in the quadratic kinetic term, as well as derivative interactions. These would 
lead to new vertices and therefore to a more complete resumation of the graphs. 



d x r = hF\r], 



(17) 
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Finally, let us make a comment on N = 2 supersymmctry. The latter can be 
obtained by dimensional reduction of N = 1 supersymmetry in 3+1 dimensions 
and thus has non-renormalization theorems. The present method can then still 
be useful for the study of the kinetic operators, since these do get renormalized. 
The evolution equation with the parameter A would be quite different though, 
due to the existence of two kinds of integrals: J d 3 xd 2 8 and J d 3 xd 2 6d 2 9. Indeed, 
the functional derivatives have to be taken with respect to the full superspace 
dependence, such that [7] 

SQ(x,e,e) _n 2 6 2 {e _ d , )S 2 { Q_- § > )§ 3 {x _ x , ) (18) 



5Q(x>,6',6) 

when Q is a chiral superfield. As a consequence, the final evolution equation is 
more involved and is planned for a future work. 

Acknowledgments 

I would like to thank the Physics Department of King's College for its support. 

Appendix A: Evolution equation for the scalar 
field 

The Euclidean action, functional of the field </>, is 

S[4>]= Jd 3 x{-\^+±rnl4> 2 + ^y (19) 

and the connected graphs generator functional, function of the source j, is 
defined by W[j] = — In Z\j] where 

2\j] = J m exp - J d 3 x. 7 (x)0(x)| . (20) 

The functional derivative of W defines the expectation value field <j): 

5W 1 

= -<<f>(x)>=4>{x), (21) 



Sj(x) Z 

where 



< 4>(x) >= J V[^{x) exp - J d? y] {y)4>{y)} . (22) 



We also have 
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4>(x)<t>(y) - \ < 4>(x)4>(y) > ■ (23) 



Inverting the relation (21) which gives <j)(x) as a function of j(x), we define the 
Legendre transform T (functional of </>) of W by 

T[<P} = W\j] - J d 3 xj(x)cj>(x). (24) 
From this definition we extract the following functional derivatives: 



' )] : -J(x) (25) 



5cj>(x) 

6 2 T ( S 2 W Y 1 



54>{x)84>(y) \Sj(x)Sj{y)J 
The evolution of W with the parameter A is, according to (23), 

d\W = -gj J d 3 x < 4> 2 (x) > 

m« f ,,, . m 2 , f ,q S 2 W ,„„. 

= -fj dx ^-^J dx WWY (26) 

To compute the evolution of T with A, we have to keep in mind that its inde- 
pendent variables are 4> and A, such that 

d x Y = 8 X W + J d 3 x-^d x j(x) - J d 3 xd x j(x)0(x) = d x W. (27) 

Combining these different results, we finally obtain the exact evolution equation 
for the proper graphs generator functional T: 



rn. 



where we used the fact that (f>(—p) = 4>{p) since 4>{x) is real. Note that, in 
Eq.(28), the integration of the operator (<5 2 r) _1 has to be understood as 



f d 3 p d 3 q ( 5 2 T \"\o ^3,3 / 



(29) 
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In the approximation where we assume that 

r = + + (30) 

it is enough to consider a constant configuration O of 4>{x), which leads to 



24^° 



( ^ V 1 = * 3 <P + g) (31) 
WiPmq) J ^ ~ p 2 + Am 2 + §^ 2 ' ^ > 

where V = <5 3 (0) is the space-time volume. The evolution equations for m 2 and 
e are obtained after expanding the right-hand side of Eq.(31) in powers of (f>o, 
which leads to: 



r d 3 P f s 2 r y 1 

J (2tt) 3 l^(p)%)J^ 



= V I constant — — ( / —z - , r „ NO ) (i 2 . 



d 3 p 



2(2tt) 3 \J (p 2 + Am 2 ) 2 / V " 
/ r d 3 p 



4(2tt) 3 V7 (P 2 + Am 2 ) 3 > Vo 



=2 



= V { C ° nStant - !^aW° + 128,A 3 / 2 m 3 ^ + - j ' (32) 

such that the identification of the coefficients of </>q and 0q m the evolution 
equation (28) leads to 

d x (Xm 2 ) = m 2 Q (l 



^ = 3^V' (33) 



We finally compute the one-loop renormalization of the potential of this 
model, using the usual Feynman rules, and starting with the action for A = 1: 

5a=i = / d 3 x + ^0 2 + |0 4 } . (34) 

The one-loop correction to the mass is given by the tadpole diagram. Taking 
into account the symmetry factor (1/2), we have 
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2 2 * e o f d 3 i 1 

2 y (2tt) 3 9 2 + to 2 , 
= m0 -R_ + »_ + O^ x j (35) 

where A is the cut-off. This last result coincides with Eq.(6) if we make the 
identification M 2 — > fte A/47r 2 . 

The one-loop correction to the charge is given by, at the limit of zero incom- 
ing momenta and taking into account the symmetry factor (1/2) as well as the 
3 different possibilities for the incoming momenta: 

. , , , ; <!■'■</ 1 



% 2 [lG) J (2^)3(^+^)2 

i 

no' 

which is consistent with the result (6) 
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Appendix B: Evolution equation for the Wess- 
Zumino model 

The reader can find a complete presentation of superspace in 2+1 dimensions 
in [7, 8] and we give here only the basic properties necessary for the present 
derivation. The conventions are those given in [7]. 

The classical action we are interested in, functional of the (real) N = 1 
superfield Q, is given by 

S[Q}= J d 5 z!^QD 2 Q + ^m Q 2 + ^Q 3 Y (37) 
where d 5 z = d 3 xd 2 9 and the spinorial derivative is 

Ot>a 

where the gamma matrices are given by 7 = a 2 ,^ 1 — ia 1 ,^ 2 = icr 3 , and 
a 1 , a 2 , cr 3 are the Pauli matrices. 

The construction of the expectation value superfield Q follows the one given 
in the Appendix A: if W is the connected graphs generator functional, function 
of the supersource P, we have 

dW 
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Following the steps of Appendix A, it is straightforward to see that the evolution 
equation of the proper graphs generator functional T is 



mo f d 3 p 2 J _ / S 2 T 



SQ(p,9)5Q(p,9) 



where we used Q(—p, 0) = Q(p, 9) since Q(z) is real. In Eq.(40), the integration 
of the operator (<5r) _1 has to be understood as 

With the following functional dependence of T: 

T[Q] = J d 5 z ^QD 2 Q + \mQ 2 + |q 3 | , (42) 



we obtain 

6 2 T 



(Am + D 2 e )(2ir) 5 (p + q)S (6 — 9') 



8Q(p,9)5Q(q,9>) ' ^ ' ^ 

+gQ{p + q,6)5 2 (6-6') 1 (43) 

where D^ g = d a — p^f^O) . We now take the inverse of (S 2 T), using the 
expansion 

(A + B)- 1 = A- 1 - A^BA- 1 +A- 1 BA- 1 BA- 1 

-A^BA^BA^BA- 1 + ... (44) 

where 

A = (Xm + D 2 pe )(2n) 3 S 3 (p + q)6 2 (9-e') 

B = gQ(p + q,9)S 2 (9-9'), (45) 

such that A is diagonal in Fourier space. We will also use the following properties 
[7]: 
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S 2 (0 1 -9 2 )S 2 (0 2 -9 1 ) = O 
S 2 (e 1 -9 2 )D« f)i 6 2 (9 2 -e 1 ) = 
2//) n \ r>2 x2/n a \ _ &l 



5\9 1 -9 2 )D z p9i 8 z {9 2 -9 1 ) = 5 



(Die? = -v\ (46) 



such that 



Am — D 2 

* ~ a 7 (2n) 3 8 3 (p + q)6 2 (9-9'). (47) 



p 2 + (Am) 5 



In the expansion (44), the constant and linear terms in B lead to constants 
after taking the trace (41). The trace of the quadratic term is 

» 2 / < 48 » 

xQ(-p + Mi)Q(-fc + ;P, 02) 

Am - D 2 H Am - £)L „ , 

X . P > 2 A - ^1 T2-TT % 5 2 (0l - 2 ) 

p z + (Am) 2 k z + (Am) 2 

Am - D 2 a „ 

x 2— n % 6 ^-9). 

p A + (Amy 

Using the properties (46), we obtain then for the quadratic term 



J (2*0 



; , <Pp__£k_f 1 

r)3 (2^)3 [p 2 + (Am) 2 ] 2 [(fc - p) 2 + (Am) 2 l 



x (3AmQ(fc, 9)D 2 _ p Q{-k, 9) + [p 2 - 3(Am) 2 ]Q(fc, 9)Q(-k, 9)) . 

In the previous integral, both terms give a kinetic contribution, when k ^ 
0. In the present approximation (42), we neglect these kinetic terms and the 
remaining contribution is obtained for k = in the second integrand. This gives 
the quantum corrections to the mass, obtained from Eq.(10). 

The trace of the cubic term gives 



-° 3 J < 50 » 

xQ(-p + q, 6 1 )Q{-q + k, 9 2 )Q(-k + p, 3 ) 
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Am — D'ic, „ Am — D 2 ft „ 

p A + (Am) 2 q z + (Am) 2 

x Am-i^ Am-!)* 

fc 2 + (Am) 2 p z + (Am) 2 

Using again the properties (46), we obtain 

5 3 /■ d 3 p d 3 q d 3 k d 2 9 



(2tt) 9 7 [p 2 + (Am) 2 ] 2 [(p + q) 2 + (Am) 2 ][(p - k) 2 + (Am) 2 ] 

x (4(Am) 3 Q(fc, 0)Q(g, 0)Q(-fc - q, 9) 
+ 6(\m) 2 Q(q, 6)Q{-k - q, 9)D 2 _ p Q(k, 9) 
+ 4AmQ(-fc - q, 9)D 2 _ p _ q [Q(q, 9)D 2 _ p Q(k, 9)] 
+p 2 Q(q, 9)Q(k, 9)D 2 _ p Q(-k - q, 9)) (51) 

The evolution of the coupling constant g is obtained in the limit of zero incom- 
ing momenta, since in the approximation (42) the derivative interactions are 
neglected. Considering the fact that 

D 2 _ p _ q D 2 _ p = -p 2 + D 2 _ q D 2 _ p - 2pq9 2 D 2 _ p , (52) 
the non-derivative interaction terms coming from Eq.(51) are then 

/d 3 a d 3 k 
Wyj2nr d20 3(M)QM)Q(-fc-<7,0) (53) 

d 3 p 4(Am)[(Am) 2 - p 2 ] 

(2tt) 3 [p 2 + (Am) 2 ] 2 [(p + q) 2 + (Am) 2 ] [(p - k) 2 + (Am) 2 ] ' 

for k = q = in the integrand. The identification of both sides of the evolution 
equation (40) gives then Eq.(15). 



We finally compute the one-loop renormalization of this model, using the 
Feynman graphs technique, and starting with the bare action for A = 1: 



(54) 



The one-loop correction to the quadratic term in Q is given by the self-energy 
graph 

r2[0] = gl j ( t]^ty d2dd29 ' Q{ - p > d)Qip >^ 
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x mo_D^^ _ m -D q+pe 
Using the properties (46), we obtain 

r 2[Q] = .9o / ff^ff^ ^ f 2^2/2x^^21 ( 56 ) 
J (2tt) 3 (2tt) 3 (mg + <7 2 )[m 2 + (9+p) 2 ] 

x {-2(? 2 (6> - 9')m Q(-p, 6)Q(p, 9) 

+ S 2 (9 9')Q(p,6>)D 2 g0 lQ(-p,6)D 2 q+p0 ,5 2 (9> 9)}} 



2 /■ d 3 p rf 3 g j2a Q(p,9)D 2 e Q(-p,0)-2m o Q(p,6)Q(-p,9) 



r) 3 (m^ + g 2 )K + ( 9 + p) 2 



The correction to the mass is obtained for p = in the non-derivative inte- 
grand, such that the one-loop mass is 



1 ) / d 3 q 2m c 

m y ' — mo ' 



R 2 M • 



) 3 (m 2 , + g 2 ) 2 

= mo-Ti^-, (57) 

47T 

and the constant c appearing in Eq.(13) is, at one loop, c = — hg 2 / (4irm ) . 
The one-loop correction to the cubic term in Q is given by 

TslQ] - 9 3 of^^^0 d We 2 (58) 
xQ(p,0)Q(q,e 1 )Q(-p-q,6 3 ) 

x i^k^ _ ^) m 2 °:^^ 2 ^ - ^ 

rriQ + k 2 m^ + (k+ P ) 2 



m 2 + (k + v + a) 2 



(k+p + q) 



The properties (46) give then 

r rni f rf 3 P d\ d 3 k d 2 8 

3M ~ (27r) 9 7 (m 2 + k 2 )[m 2 + (k+p) 2 ][m 2 + (k+p + q) 2 ] 
x {3m Q(-p - q, 6)D 2 k+pe [Q( Pl 6)Q(q, 9)] 

- 3m 2 Q(p,9)Q(q,9)Q(-p - q,9) 

- D 2 k+p9 [Q(q, 9)Q(-p - q, 9)D 2 g Q{P, 0)]} 
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(k 2 -3mt)Q(p,6)Q(q,9)Q(-p-q,6) 



{m 2 + k 2 )[ml + {k + p) 2 ] [m 2 + (k + p + q) 2 ] 



+ derivative terms 



(59) 



The correction to the coupling constant is obtained in the limit p, q — ► in the 
non-derivative integrand, i.e. it is given by the integral 



such that the one-loop correction to the coupling constant vanishes, in accor- 
dance with the the result (16). 
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